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Abstract 

We compute the photon number distribution, the Q(a) distribution function and 
the wave functions in the momentum and position representation for a single 
mode squeezed number state using generating functions which allow to obtain 
any matrix element in the squeezed number state representation from the matrix 
elements in the squeezed coherent state representation. For highly squeezed 
number states we discuss the previously unnoted oscillations which appear in 
the Q(ct) function. We also note that these oscillations can be related to the 
photon-number distribution oscillations and to the momentum representation 
of the wave function. 
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1 Introduction 



The study and experimental detection of squeezed states jl], ^, |^, || and other 
non classical states of optical systems, is an interesting issue, from both the 
fundamental and the technological point of view. On one hand, disentangling the 
properties of intrinsically quantum states of light enhances our understanding 
of the behavior and interactions of photons, and improves our knowledge of 
Quantum Electrodynamics. On the other hand, the use of light states with 
reduced quantum fluctuations in one of the conjugate quadratures is at the 
heart of attractive proposals for the detection of very weak signals H p). 

In this article we are concerned with the study of the squeezed number states 
and their phase space representations . These states emerged in the original 
analysis of Yuen |l) of the Two Photon Coherent States and have received since 
then a moderate amount of attention in the literature [Q, [|, ||, |l0|; [llj [l^] . They 
appear promissory for discussing the departure from the semiclassical regime in 
the context of quantum optics and for this reason were chosen for this study, 
but they also may have a practical interest on their own ]T^ | . 

The one mode squeezed number states are obtained by applying the squeeze 
operator to the Fock states. As we discuss below, their wave functions in the 
position and momentum like quadratures are obtained by literally compressing 
or stretching the corresponding wave functions of the Fock states. An impor- 
tant aspect which holds our attention through this work is that they present 
oscillations in the photon number distribution as have been reported in Refs. 
[0 Pp tT These oscillations were first observed in the photon number distri- 
butions of squeezed states |l3| which are equivalent to the Two-Photon Co- 
herent States, and present several interesting properties. In a series of papers, 
[ |i"3[ |l4| , [l5| |l6| , [TtJ , these oscillations were explained as the result of the interfer- 
ence of different contributions to probability amplitudes. These contributions 
were interpreted as the overlapping of the regions assigned to each quantum 
state on the phase space. For the squeezed number states oscillations, a closely 
related discussion was sketched in [0, ||. This treatment which is physically 
very appealing is based on the semiclassical representation of states as bounded 
areas in phase space obtained with the help of Wigner distributions (l^] . In 
Ref. for the case of the squeezed states this area was fixed by weighting the 
Bohr-Sommerfeld rings with the Wigner-Cohen [po| distribution function. The 
same prescription was followed in Refs. |t], |§|, ^) for the squeezed number states. 
In Ref ]l8| , using the WKB method, a more general framework for the use of 
this picture was established. 

In this paper we study the properties of the squeezed number states with the 
help of various representations. We compute explicitly the probability ampli- 
tudes and distributions of the squeezed number states in the Fock, position and 
momentum representations. We also compute the Q(a) distribution which has 
the advantage of being a true probability distribution which can be interpreted 
as the probability of a simultaneous measurement of position and momentum 
within some region of the phase space. For high squeezing, the Q(a) function of 
a squeezed number state presents intrinsic oscillations in the space of the com- 
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plex labels a. These oscillations were not discussed previously in the literature 
and can be correlated with those in the momentum and Fock representations. 

Our computations are done using a generating function technique which can 
be applied quite generally to compute expectation values of arbitrary observ- 
ables in different representations. 

2 Two photon coherent states and squeezed num- 
ber states 

We consider here a single mode of the radiation field described in terms of 
the creation and annihilation operators a', a, the number operator N a = a^a, 
and the momentum-like and position- like quadratures p = ~^( a ~ an d q = 
■^=(a+(r). Number (Fock) states \n) and coherent states \(3) are the eigenstates 

of N a and a with eigenvalues n and (3 respectively. Coherent states may also 
be characterized as displaced vacuum states and are minimal uncertainty states 
satisfying Aq 2 — Ap 2 = i. We are working in units such that h — 2ir. In the 
number state representation they are written as: 



= e -mv^^\n). (1) 



Squeezed states and squeezed number states may be obtained directly from \(3) 
and \n) by the application of the squeeze operator. This operator depends on 
a complex parameter £ = re~ l< ^. For simplicity we restrict our computation to 
the case of real £. We define S{r) as the squeeze operator, and therefore 

\p,r) = S(r)\P) (2) 

|m, r) — S(r)\m) (3) 

where S(r) is given by 

S(r)=expQr(at) 2 -Ira 2 ). (4) 

It is also useful to define a transformed annihilation operator b as: 

b = SaS^ = cosh(r)a + sinh(r)a 1 ' (5) 

and then, one has 

b\/3,r) = p\p,r) (6) 

and 

N})\m,r) = b^b\m^r) = m\rn,r) . (7) 

Given the linear nature of the transformation on a (||) , one readily computes 
the coherent states and Fock states amplitudes for the squeezed states [|l|, 
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(a\0,r) 



cosh 1 ^ 2 (r) 



exp 



tanh(r) 



tanh(r) 2 



cosh(r) 



Q*/3 



(8) 



tanh"/ 2 (r) f /3 



(2"n! cosh(r))V2 « ^ ( 2 s i n h(r) cosh(r)) 1 /2 

xexp(-i|/f + itanh(r)/3 2 ) (9) 

where H ra is the nth order Hermite polynomial. 

For /3 e M, the wave functions in 5 and p are then given by: 

MM = (2 7 rA 2 <Z )- 1 /4 exp (-^|^) (10) 

(p|A = (^) V4 cxp (-p 2 A 2 q - ipq ) , (11) 
with g = v / 2e r /3, A 2 q = e - 2r /2 and A 2 p = e 2 ''/2. 



3 Representations of the squeezed number states 



To discuss the properties of the squeezed number states, we observe that for 
any operator R and any state \4>) , the completeness of the Fock states and Eq. 
(Q) allow one to write, 



X>|AS(r)|» 



= e -l/3|V2 



E 



P m (<ip\R\m,r) 
(ml) 1 / 2 



(12) 



This provides a generating function for the matrix element (ip\R\m, r) which 
can then be obtained as, 



\R\m, r) 



1 



(ml) 1 / 2 



d" 



dp 



-{e^ 2 ' 2 ^\R\P,r)) 



(13) 



/3=0 



For example, if R = U(t) is the evolution operator, and = \k) with 
{|fc)} being any complete basis of the radiation states, one can compute in this 
form the time dependent probability amplitudes in the {|fc)} representation. For 
t = we get the probability amplitudes of stationary squeezed number states. 

Generalizing the same idea, we also note that 



(a g , r\R\/3 g , r) = e 



-1/2M 2 -1/2|/3| ; 



EE 



(nl) 1 / 2 (ml) 1 / 2 



(n,r\R\m,r) (14) 
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so that we can write, 



(n, r\R\m, r) 



1 



gn gr, 



(n!m!) 1 /2 da* n df3 r - 



3 l/2|a| 2 1/2|/3| 2 



'■' " « 1 (a s ,r|ii|/3 ff ,r) 

a*=0./3=0 

(15) 

With this observation we can generate the probability distributions for the 
squeezed number states in the Fock, position, momentum, and coherent states 
representations. Taking — \n) and R = 1, the amplitude for the squeezed 
number states in the Fock states basis is, 



(n|ro, r) 



1 



(ml) 1 / 2 



d" 



(tanh(r)/2)"/ 2 
\Jn\m\ cosh(r) 

Next we note that, 



9/3' 



-H„ 







-(e^ 2 / 2 (n|/3 g ,r)) 



0=0 



y/ 2 sinh(r) cosh(r) 



f3=o 



H W 



(2 sinh(r) cosh(r)) 1 / 2 



0=0 



2'(_1)«-V2 n !( 2s inh(r) cosh(r))- z /2 



and 







(e 2 ) W |/3=o 

Using the Cauchy formula 

{f(x)g(x)}W = ^ 

we get, 
(n|m, r) 



/ tanh(r) \ 

V 2 J 



i/2 



I odd 
i even 



A//! 



fc=0 



k\{m-k)\ 



min(m,n) 

cosh(r) 2 fc 



(win!) 1 / 2 ; .sinh(r) . (n+m-a>) (— 1) : 



-)- 



2 J'\ 2 

0,2,4,6... n, m even 
1, 3, 5, 7... ?i, to odd 



(16) 



(17) 



(18) 



(19) 



(20) 



This result is in agreement with the one reported in Ref . |7| , which was ob- 
tained by the application of a normally ordered squeeze operator to a Fock state. 
The most salient features of this photon distribution are that it oscillates and 
that only the photon number states of the same parity of m are represented in 
the expansion for \m,r). The latter reflects of course the quadratic dependence 
of S(r) in a and a'. The oscillations of the distribution P n , m = | (n\m, r) | 2 for 
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m > 1 and r sufficiently large have a fixed number of maxima; m/2 + 1 when 
m is even and (to + l)/2 when m is odd. These characteristics can be related 
to the structure of the Q(a) function, as we discuss below. 

For the squeezed number states with to = (the squeezed vacuum) and 
m = 1 the sum in equation (20) is just a single term: 

, , , . , 9 n\ tanh"(r) 

^°' r)l = ((|)!)^ cosh(r) (21) 
n = 0,2,4,6, .... 

i/ i- m9 Til tanh n_1 (r) ,„ . 

■< n|1 ' r)l = a^)!)^ cosh 3 (r ) ■ ^ 

n = 1, 3, 5, 7, .... 

As an example, we show in Figure ([l]) the photon number distribution for a 
squeezed number state with m — 7 and r = 1.4. The distribution |(n|7, 1.40) | 2 
has 4 maxima in n located at 1,11,37 and 89 photons respectively. 

The support of the distribution, and the oscillations, become wider for 
greater r. For the squeezed number state with to = 7, this is illustrated by 
Figure (|^). In this figure, we plot the positions of n for the last three maxima 
as functions of r. These values shift to the right as r grows. It is possible 
to observe a tendency for each of the maxima, which can be approximated by 
exponential functions . 

Let us study the other representations of squeezed number states. Using 
(0) and (H ) we have 

am ( ( -{g-V2^fi) 2 \ 

(g|m,r)=m!- 1 / 2 (2 7 r)- 1 / 4 (2- 1 /2 e - ) -t/2^_n — ) ^ 

"^(23) 

And using again the Cauchy formula we have 

6 J X — ^^~ r — ) I 2 \ = V e\ — — ) — e 13 I 2 

d(3 m \ J k)\k\df3 k df3 m - k 

Evaluating at (3 = we obtain 

(,|,„, r) . .-'V^),-VE ^ (0) ,25) 

0, 2, 4, 6... to even 
1, 3, 5, 1 ... to odd 

In order to simplify the expression above we use the following identity of the 
Hermite polynomials, 

1 m 2 k/2 ( x \ 

^! Hm(x) = X- fc!((TO-fc)/2)! Hfc \WH) (26) 



(24) 
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0. 2. 4. 6... m even 

1, 3, 5, 7... m odd 

Finally, the amplitude (^5|) becomes: 

(q\m, r) = ^/V^^-A^) 2~ m ' 2 m\~ 1 l 2 ^ m (-^) (27) 

Thus, the amplitude is that of a Fock state, depending on a squeezed quadrature 
variable e r q. Moreover, the amplitude of a Fock state is given by the limit of 
( |7|) when r — > 0. 

The representation of the squeezed number states in terms of the other 
quadrature follows along the same lines. We get 

0, 2, 4, 6... m even 

1, 3, 5, 1 ... m odd 

and 

<PK r ) = JyC/A '^^ 2- m/2 HrH m ( e -» (29) 

The amplitude depends on e~~ r p and, for r > 0, |(p|to, r)| 2 has a support that 
is e r wider than the support of a Fock state. The distribution has m+1 maxima 
and m minima (which are also zeroes). In Figure ([?]) we show the momentum 
probability distribution for the m = 7 squeezed number state with r — 1.4. The 
oscillations in this distribution come from the Hermite polynomial. For large 
squeezing, these oscillations reappear in the Q(a) distribution function. 

Finally, let us turn to the Q(a) function, which is given by the diagonal 
elements of the density matrix p associated with the state in the coherent states 
basis. 

For a squeezed number state one has, 

Q(a) = -(a\p\a) = -\(a\m,0\ 2 - (30) 

The function Q{a) gives the probability of being in a minimum dispersion 
state around an average position and momentum proportional to Re(a) and 
Im(a) pl| . Following the previous procedure and using (|8|), we have 

Hm , r> = ( \ y /2 e -l/2|a|'- ta nh(r)a-V2 ^_ e ^P' + ^,P 

v 1 ' ' V m!cosh ( r ) / d P m 

(31) 

which can be reduced to 

(a\m,r) = m\ 1 ' 2 cosh(r)- 1 / 2 e - 1 / 2 l«l 2 -t^(r)«« 2 /2 



G 



[W2] n -r 

SI 

(m-2p)\p\ 

m/2 to even 
^2^!- to odd 



[m/2] = 

The Q{a) function of a squeezed number state is, therefore: 
Q(a) = -m\cosh(r)- 1 e- lal2 - t!inhir)(a ' 2+a2)/2 



7T 



^ 2-?sintr>) cosher) p 



p=0 



(rn — 2p) !p 
[to/2] 



(33) 



to/2 to even 
^ to odd 

The shape of Q(a) depends strongly on the squeezing parameter r, as illus- 
trated by Figures (@), (|) and (|), (0); which show the Eq. © for TO = 7. 
r = 0.50 and r = 1.40. 

For small squeezing, the Q(a) function is appreciably nonzero over an ellip- 
tical ring, with an eccentricity e r . When r = this ring is circular, the Q{a) 
function being that of a Fock state. For greater squeezing , the Q{a) function 
shows prominent oscillations with a number of maxima which for the to squeezed 
number state stabilize at to + 1. As can be shown by direct numerical computa- 
tion, the transition occurs at r ~ i ln(m) which corresponds to y/me~ r = 1. It 
is interesting to note that in the phase space analysis this behavior appears at 
the point where the inner boundary of the deformed ring touch itself and in this 
sense is also a consequence of interference in the phase space. Moreover it can 
be checked by direct numerical computation that this effect is registered only 
by the Q function, and not by other interesting distributions like the Wigner or 
Wigner-Cohen functions. 

The oscillations for r > | ln(m) occur along the axis of Im(a) and are corre- 
lated to the oscillations of the momentum wave function. To show this explicitly 
observe that for high squeezing, the Q(a) function evaluated on the Im(a) axis 
is proportional to the square of the momentum wave function: 

Q(a = iV2p)~\ (p\ to, r) | 2 (34) 

For illustration, consider the to — 7 squeezed number state with r = 1.4. The 
maxima in its momentum distribution (see Figure (||)) correspond to those of 
the Q(a) function (figure ^|) in the Im(a) axis. 

In the same way, for large squeezing, the Q(a) oscillations are also corre- 
lated to the oscillations of the photon number distribution |(n|TO,r)| 2 . Each 
maximum a max = i \a max \ (\a max \ > 1) of Q is associated to a maximum n max 
of |(n|TO,r)| 2 , in such away that, 

Tlmax — \(%max\ • 
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In this form we can use the Q function to obtain directly from the phase 
space a good picture of P m distribution of the state. 

Finally, observe that the Q distribution of a squeezed number state is given 
by, 

Q(a) = -\< a\S\m > | 2 = -| < m\S*\a > | 2 . 

and henceforth, the Q(a) function of a squeezed number state \m,r) is pro- 
portional to the photon number amplitude of a two photon coherent state 
\a,r) = S^\a >. The oscillations still can be explained using the hypothesis 
of interference in phase space. To illustrate this point consider Figure (0) which 
shows the areas representing a Fock state and a pair of two photon coherent 
states ( |ck,r)) with respectively real or imaginary a . For a imaginary the 
overlapping area doubles and produces the interference effect in the Q function. 
The photon statistical value of a two photon coherent state is computed using 
the method of overlapping areas in the phase space j2^] . The result is, 

P m = | < m\S*\a = iy > | 2 = \\fAme i<l> + vCW"^| 2 (35) 

with 

e -2(m+l/2-y 2 /e 2r )/e 2r 
An = , (36) 

y/2n(m + l/2-y 2 /e 2r ) 

and 



(j) = (m + 1/2) arctan{^(m + 1/2 - y 2 / e 2r / (y 2 / e 2r )} 

-xy/(m + 1/2 - y 2 /e 2r /e 2r ) - tt/4 (37) 

This match with the exact computation presented above as is shown in Figure 
(U) where we present respectively the Q(a — iy) function for the m = 7, r = 1.4 
squeezed number state obtained using this approximation and the previous exact 
calculation. For sufficiently high values of \a\ the WKB approximation on which 
is based the phase space computation no longer holds as evidenced in Figure 



4 Conclusion 

We have computed the photon number distribution, the momentum like and po- 
sition like wave functions, and the Q(a) function for squeezed number states and 
we have shown that each of them has characteristic oscillations which depend 
on the squeeze parameter r. For highly squeezed number states we observed 
that the oscillations in the different probability distributions are fixed in num- 
ber and are in close correspondence between them. The Q(a) function in this 
case present a richer structure that for number states. The oscillations were also 
discussed within the semiclassical approximation using the phase space methods 
based in the WKB functions (lj, ||. 
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It should be also mentioned that the generating function formalism used 
throughout this work has been shown to be a straightforward and very valuable 
procedure for calculating probability amplitudes and matrix elements for any 
base obtained transforming number states by any arbitrarily chosen unitary 
operator. 
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Figure 1: Photon number distribution for the squeezed number states with 
to = 7 and r = 1.40, \(n\7, 1.40)| 2 . 



11 



100 



80 



60 
n 



40 



20 



0.9 1 1.1 1.2 1.3 1.4 1.5 

r 



Figure 2 : Position n of the maxima for the photon number distribution | (n 1 7, r) \ 2 
of a squeezed number states, versus r 
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Figure 3: Momentum probability distribution | (p\m, r) | 2 for the m = 7 squeezed 
number state with r = 1.4. 
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Figure 4: Q(a) function of the m — 7 squeezed number state with r = 0.50 
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Figure 5: Contour graphic of the Q{a) function of the m = 7 squeezed number 
state with r = 0.50 
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Figure 6: Q{a) function of the m — 7 squeezed number state with r = 1.40 
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Figure 7: Contour graphic of the Q{a) function of the m = 7 squeezed number 
state with r = 1.40 
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Figure 8: Phase space areas for: a number state \m > (circular ring), a two 
photon coherent states S^\a > with real a ( right bottom ellipse ) and a two 
photon coherent states with imaginary a (center top ellipse) 
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Figure 9: Q(a = iy) function for the m = 7, r = 1.4 squeezed number state: 
The two areas overlapping approximation ( left) and the exact calculation ( 
right). 



18 



